The architecture of redox microdomains: Cascading gradients and peroxiredoxins' redox-oligomeric coupling integrate redox signaling and antioxidant protection.

In the cytosol of human cells under low oxidative loads, hydrogen peroxide is confined to microdomains around its supply sites, due to its fast consumption by peroxiredoxins. So are the sulfenic and disulfide forms of the 2-Cys peroxiredoxins, according to a previous theoretical analysis [Travasso et al., Redox Biology 15 (2017) 297]. Here, an extended reaction-diffusion model that for the first time considers the differential properties of human peroxiredoxins 1 and 2 and the thioredoxin redox cycle predicts important new aspects of the dynamics of redox microdomains. The peroxiredoxin 1 sulfenates and disulfides are more localized than the corresponding peroxiredoxin 2 forms, due to the former peroxiredoxin's faster resolution step. The thioredoxin disulfides are also localized. As the H2O2 supply rate (vsup) approaches and then surpasses the maximal rate of the thioredoxin/thioredoxin reductase system (V), these concentration gradients become shallower, and then vanish. At low vsup the peroxiredoxin concentration determines the H2O2 concentrations and gradient length scale, but as vsup approaches V, the thioredoxin reductase activity gains influence. A differential mobility of peroxiredoxin disulfide dimers vs. reduced decamers enhances the redox polarity of the cytosol: as vsup approaches V, reduced decamers are preferentially retained far from H2O2 sources, attenuating the local H2O2 buildup. Substantial total protein concentration gradients of both peroxiredoxins emerge under these conditions, and the concentration of reduced peroxiredoxin 1 far from the H2O2 sources even increases with vsup. Altogether, the properties of 2-Cys peroxiredoxins and thioredoxin are such that localized H2O2 supply induces a redox and functional polarization between source-proximal regions (redox microdomains) that facilitate peroxiredoxin-mediated signaling and distal regions that maximize antioxidant protection.

Pillars of theoretical biology: "Biochemical systems analysis, I, II and III".

Michael Savageau's Biochemical Systems Analysis I, II, IIIpapers, published in volumes 25 and 26 of the journal,kickstarted a research programme that originated many of the core concepts and tools of Systems Biology. This article briefly summarizes these papers anddiscusses the most relevant developments in Biochemical Systems Theory since their publication.

The mammalian peroxisomal membrane is permeable to both GSH and GSSG - Implications for intraperoxisomal redox homeostasis.

Despite the large amounts of H2O2 generated in mammalian peroxisomes, cysteine residues of intraperoxisomal proteins are maintained in a reduced state. The biochemistry behind this phenomenon remains unexplored, and simple questions such as "is the peroxisomal membrane permeable to glutathione?" or "is there a thiol-disulfide oxidoreductase in the organelle matrix?" still have no answer. We used a cell-free in vitro system to equip rat liver peroxisomes with a glutathione redox sensor. The organelles were then incubated with glutathione solutions of different redox potentials and the oxidation/reduction kinetics of the redox sensor was monitored. The data suggest that the mammalian peroxisomal membrane is promptly permeable to both reduced and oxidized glutathione. No evidence for the presence of a robust thiol-disulfide oxidoreductase in the peroxisomal matrix could be found. Also, prolonged incubation of organelle suspensions with glutaredoxin 1 did not result in the internalization of the enzyme. To explore a potential role of glutathione in intraperoxisomal redox homeostasis we performed kinetic simulations. The results suggest that even in the absence of a glutaredoxin, glutathione is more important in protecting cysteine residues of matrix proteins from oxidation by H2O2 than peroxisomal catalase itself.

Ligand's Partition to the Lipid Bilayer Should Be Accounted for When Estimating Their Affinity to Proteins.

Ligand-protein interactions are usually studied in complex media that also contain lipids. This is particularly relevant for membrane proteins that are always associated with lipid bilayers, but also for water-soluble proteins studied in in vivo conditions. This work addresses the following two questions: (i) How does the neglect of the lipid bilayer influence the apparent ligand-protein affinity? (ii) How can the intrinsic ligand-protein affinity be obtained? Here we present a framework to quantitatively characterize ligand-protein interactions in complex media for proteins with a single binding site. The apparent affinity obtained when following some often-used approximations is also explored, to establish these approximations' validity limits and to allow the estimation of the true affinities from data reported in literature. It is found that an increase in the ligand lipophilicity or in the volume of the lipid bilayer always leads to a decrease in the apparent ligand-protein affinity, both for water-soluble and for membrane proteins. The only exceptions are very polar ligands (excluded from the lipid bilayer) and ligands whose binding affinity to the protein increases supralinearly with ligand lipophilicity. Finally, this work discusses which are the most relevant parameters to consider when exploring the specificity of membrane proteins.

How abundant are superoxide and hydrogen peroxide in the vasculature lumen, how far can they reach?

Paracrine superoxide (O2•-) and hydrogen peroxide (H2O2) signaling critically depends on these substances' concentrations, half-lives and transport ranges in extracellular media. Here we estimated these parameters for the lumen of human capillaries, arterioles and arteries using reaction-diffusion-advection models. These models considered O2•- and H2O2 production by endothelial cells and uptake by erythrocytes and endothelial cells, O2•- dismutation, O2•- and H2O2 diffusion and advection by the blood flow. Results show that in this environment O2•- and H2O2 have half-lives <60. ms and <40. ms, respectively, the former determined by the plasma SOD3 activity, the latter by clearance by endothelial cells and erythrocytes. H2O2 concentrations do not exceed the 10 nM scale. Maximal O2•- concentrations near vessel walls exceed H2O2's several-fold when the latter results solely from O2•- dismutation. Cytosolic dismutation of inflowing O2•- may thus significantly contribute to H2O2 delivery to cells. O2•- concentrations near vessel walls decay to 50% of maximum 12 µm downstream from O2•- production sites. H2O2 concentrations in capillaries decay to 50% of maximum 22 µm (6.0 µm) downstream from O2•- (H2O2) production sites. Near arterioles' (arteries') walls, they decay by 50% within 6.0 µm (4. µm) of H2O2 production sites. However, they reach maximal values 50 µm (24 µm) downstream from O2•- production sites and decrease by 50% over 650 µm (500 µm). Arterial/olar endothelial cells might thus signal over a mm downstream through O2•--derived H2O2, though this requires nM-sensitive H2O2 transduction mechanisms.

Analysis of the Equilibrium Distribution of Ligands in Heterogeneous Media-Approaches and Pitfalls.

The equilibrium distribution of small molecules (ligands) between binding agents in heterogeneous media is an important property that determines their activity. Heterogeneous systems containing proteins and lipid membranes are particularly relevant due to their prevalence in biological systems, and their importance to ligand distribution, which, in turn, is crucial to ligand's availability and biological activity. In this work, we review several approaches and formalisms for the analysis of the equilibrium distribution of ligands in the presence of proteins, lipid membranes, or both. Special attention is given to common pitfalls in the analysis, with the establishment of the validity limits for the distinct approaches. Due to its widespread use, special attention is given to the characterization of ligand binding through the analysis of Stern-Volmer plots of protein fluorescence quenching. Systems of increasing complexity are considered, from proteins with single to multiple binding sites, from ligands interacting with proteins only to biomembranes containing lipid bilayers and membrane proteins. A new formalism is proposed, in which ligand binding is treated as a partition process, while considering the saturation of protein binding sites. This formalism is particularly useful for the characterization of interaction with membrane proteins.

Modeling Gd3+ Complexes for Molecular Dynamics Simulations: Toward a Rational Optimization of MRI Contrast Agents.

The correct parametrization of lanthanide complexes is of the utmost importance for their characterization using computational tools such as molecular dynamics simulations. This allows the optimization of their properties for a wide range of applications, including medical imaging. Here we present a systematic study to establish the best strategies for the correct parametrization of lanthanide complexes using [Gd(DOTA)]- as a reference, which is used as a contrast agent in MRI. We chose the bonded model to parametrize the lanthanide complexes, which is especially important when considering the study of the complex as a whole (e.g., for the study of the dynamics of its interaction with proteins or membranes). We followed two strategies: a so-called heuristic approach employing strategies already published by other authors and another based on the more recent MCPB.py tool. Adjustment of the Lennard-Jones parameters of the metal was required. The final topologies obtained with both strategies were able to reproduce the experimental ion to oxygen distance, vibrational frequencies, and other structural properties. We report a new strategy to adjust the Lennard-Jones parameters of the metal ion in order to capture dynamic properties such as the residence time of the capping water (τm). For the first time, the correct assessment of the τm value for Gd-based complexes was possible by recording the dissociative events over up to 10 µs all-atom simulations. The MCPB.py tool allowed the accurate parametrization of [Gd(DOTA)]- in a simpler procedure, and in this case, the dynamics of the water molecules in the outer hydration sphere was also characterized. This sphere was divided into the first hydration layer, an intermediate region, and an outer hydration layer, with a residence time of 18, 10 and 19 ps, respectively, independent of the nonbonded parameters chosen for Gd3+. The Lennard-Jones parameters of Gd3+ obtained here for [Gd(DOTA)]- may be used with similarly structured gadolinium MRI contrast agents. This allows the use of molecular dynamics simulations to characterize and optimize the contrast agent properties. The characterization of their interaction with membranes and proteins will permit the design of new targeted contrast agents with improved pharmacokinetics.

Glucose-6-Phosphate Dehydrogenase Deficiency and Neonatal Hyperbilirubinemia: Insights on Pathophysiology, Diagnosis, and Gene Variants in Disease Heterogeneity.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a prevalent condition worldwide and is caused by loss-of-function mutations in the G6PD gene. Individuals with deficiency are more susceptible to oxidative stress which leads to the classical, acute hemolytic anemia (favism). However, G6PD deficiency in newborn infants presents with an increased risk of hyperbilirubinemia, that may rapidly escalate to result in bilirubin induced neurologic dysfunction (BIND). Often with no overt signs of hemolysis, G6PD deficiency in the neonatal period appears to be different in the pathophysiology from favism. This review discusses and compares the mechanistic pathways involved in these two clinical presentations of this enzyme disorder. In contrast to the membrane disruption of red blood cells and Heinz bodies formation in favism, G6PD deficiency causing jaundice is perhaps attributed to the disruption of oxidant-antioxidant balance, impaired recycling of peroxiredoxin 2, thus affecting bilirubin clearance. Screening for G6PD deficiency and close monitoring of affected infants are important aspects in neonatal care to prevent kernicterus, a permanent and devastating neurological damage. WHO recommends screening for G6PD activity of all infants in countries with high prevalence of this deficiency. The traditional fluorescent spot test as a screening tool, although low in cost, misses a significant proportion of cases with moderate deficiency or the partially deficient, heterozygote females. Some newer and emerging laboratory tests and diagnostic methods will be discussed while developments in genomics and proteomics contribute to increasing studies that spatially profile genetic mutations within the protein structure that could predict their functional and structural effects. In this review, several known variants of G6PD are highlighted based on the location of the mutation and amino acid replacement. These could provide insights on why some variants may cause a higher degree of phenotypic severity compared to others. Further studies are needed to elucidate the predisposition of some variants toward certain clinical manifestations, particularly neonatal hyperbilirubinemia, and how some variants increase in severity when co-inherited with other blood- or bilirubin-related genetic disorders.

Calculation of Permeability Coefficients from Solute Equilibration Dynamics: An Assessment of Various Methods.

Predicting the rate at which substances permeate membrane barriers in vivo is crucial for drug development. Permeability coefficients obtained from in vitro studies are valuable for this goal. These are normally determined by following the dynamics of solute equilibration between two membrane-separated compartments. However, the correct calculation of permeability coefficients from such data is not always straightforward. To address these problems, here we develop a kinetic model for solute permeation through lipid membrane barriers that includes the two membrane leaflets as compartments in a four-compartment model. Accounting for solute association with the membrane allows assessing various methods in a wide variety of conditions. The results showed that the often-used expression Papp = ß × r/3 is inapplicable to very large or very small vesicles, to moderately or highly lipophilic solutes, or when the development of a significant pH gradient opposes the solute's flux. We establish useful relationships that overcome these limitations and allow predicting permeability in compartmentalised in vitro or in vivo systems with specific properties. Finally, from the parameters for the interaction of the solute with the membrane barrier, we defined an intrinsic permeability coefficient that facilitates quantitative comparisons between solutes.

Intra-dimer cooperativity between the active site cysteines during the oxidation of peroxiredoxin 2.

Peroxiredoxin 2 (Prdx2) and other typical 2-Cys Prdxs function as homodimers in which hydrogen peroxide oxidizes each active site cysteine to a sulfenic acid which then condenses with the resolving cysteine on the alternate chain. Previous kinetic studies have considered both sites as equally reactive. Here we have studied Prdx2 using a combination of non-reducing SDS-PAGE to separate reduced monomers and dimers with one and two disulfide bonds, and stopped flow analysis of tryptophan fluorescence, to investigate whether there is cooperativity between the sites. We have observed positive cooperativity when H2O2 is added as a bolus and oxidation of the second site occurs while the first site is present as a sulfenic acid. Modelling of this reaction showed that the second site reacts 2.2 ± 0.1 times faster. In contrast, when H2O2 was generated slowly and the first active site condensed to a disulfide before the second site reacted, no cooperativity was evident. Conversion of the sulfenic acid to the disulfide showed negative cooperativity, with modelling of the exponential rise in tryptophan fluorescence yielding a rate constant of 0.75 ± 0.08 s-1 when the alternate active site was present as a sulfenic acid and 2.29 ± 0.08-fold lower when it was a disulfide. No difference in the rate of hyperoxidation at the two sites was detected. Our findings imply that oxidation of one active site affects the conformation of the second site and influences which intermediate forms of the protein are favored under different cellular conditions.

Quantitative Assessment of Methods Used To Obtain Rate Constants from Molecular Dynamics Simulations-Translocation of Cholesterol across Lipid Bilayers.

Accurately calculating rate constants of macroscopic chemical processes from molecular dynamics simulations is a long-sought but elusive goal. The problem is particularly relevant for processes occurring in biological systems, as is the case for ligand-protein and ligand-membrane interactions. Several formalisms to determine rate constants from easily accessible free-energy profiles [Δ Go( z)] of a molecule along a coordinate of interest have been proposed. However, their applicability for molecular interactions in condensed media has not been critically evaluated or validated. This work presents such evaluation and validation and introduces improved methodology. As a case study, we have characterized quantitatively the rate of translocation of cholesterol across 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine bilayers. Translocation across lipid bilayers is the rate-limiting step in the permeation of most drugs through biomembranes. We use coarse-grained molecular dynamics simulations and different kinetic formalisms to calculate this rate constant. A self-consistent test of the applicability of various available formalisms is provided by comparing their predictions with the translocation rates obtained from actual events observed in long unrestrained simulations. To this effect, a novel procedure was used to obtain the effective rate constant, based on an analysis of time intervals between transitions among different states along the reaction coordinate. While most tested formalisms lead to results in reasonable agreement (within a factor of 5) with this effective rate constant, the most adequate one is based on the explicit relaxation frequencies from the transition state in the forward and backward directions along the reaction coordinate.

Mapping the phenotypic repertoire of the cytoplasmic 2-Cys peroxiredoxin - Thioredoxin system. 1. Understanding commonalities and differences among cell types.

The system (PTTRS) formed by typical 2-Cys peroxiredoxins (Prx), thioredoxin (Trx), Trx reductase (TrxR), and sulfiredoxin (Srx) is central in antioxidant protection and redox signaling in the cytoplasm of eukaryotic cells. Understanding how the PTTRS integrates these functions requires tracing phenotypes to molecular properties, which is non-trivial. Here we analyze this problem based on a model that captures the PTTRS' conserved features. We have mapped the conditions that generate each distinct response to H2O2 supply rates (vsup), and estimated the parameters for thirteen human cell types and for Saccharomyces cerevisiae. The resulting composition-to-phenotype map yielded the following experimentally testable predictions. The PTTRS permits many distinct responses including ultra-sensitivity and hysteresis. However, nearly all tumor cell lines showed a similar response characterized by limited Trx-S- depletion and a substantial but self-limited gradual accumulation of hyperoxidized Prx at high vsup. This similarity ensues from strong correlations between the TrxR, Srx and Prx activities over cell lines, which contribute to maintain the Prx-SS reduction capacity in slight excess over the maximal steady state Prx-SS production. In turn, in erythrocytes, hepatocytes and HepG2 cells high vsup depletes Trx-S- and oxidizes Prx mainly to Prx-SS. In all nucleated human cells the Prx-SS reduction capacity defined a threshold separating two different regimes. At sub-threshold vsup the cytoplasmic H2O2 concentration is determined by Prx, nM-range and spatially localized, whereas at supra-threshold vsup it is determined by much less active alternative sinks and µM-range throughout the cytoplasm. The yeast shows a distinct response where the Prx Tsa1 accumulates in sulfenate form at high vsup. This is mainly due to an exceptional stability of Tsa1's sulfenate. The implications of these findings for thiol redox regulation and cell physiology are discussed. All estimates were thoroughly documented and provided, together with analytical approximations for system properties, as a resource for quantitative redox biology.

Localized redox relays as a privileged mode of cytoplasmic hydrogen peroxide signaling.

Hydrogen peroxide (H2O2) is a key signaling agent. Its best characterized signaling actions in mammalian cells involve the early oxidation of thiols in cytoplasmic phosphatases, kinases and transcription factors. However, these redox targets are orders of magnitude less H2O2-reactive and abundant than cytoplasmic peroxiredoxins. How can they be oxidized in a signaling time frame? Here we investigate this question using computational reaction-diffusion models of H2O2 signaling. The results show that at H2O2 supply rates commensurate with mitogenic signaling a H2O2 concentration gradient with a length scale of a few tenths of µm is established. Even near the supply sites H2O2 concentrations are far too low to oxidize typical targets in an early mitogenic signaling time frame. Furthermore, any inhibition of the peroxiredoxin or increase in H2O2 supply able to drastically increase the local H2O2 concentration would collapse the concentration gradient and/or cause an extensive oxidation of the peroxiredoxins I and II, inconsistent with experimental observations. In turn, the local concentrations of peroxiredoxin sulfenate and disulfide forms exceed those of H2O2 by several orders of magnitude. Redox targets reacting with these forms at rate constants much lower than that for, say, thioredoxin could be oxidized within seconds. Moreover, the spatial distribution of the concentrations of these peroxiredoxin forms allows them to reach targets within 1 µm from the H2O2 sites while maintaining signaling localized. The recruitment of peroxiredoxins to specific sites such as caveolae can dramatically increase the local concentrations of the sulfenic and disulfide forms, thus further helping these species to outcompete H2O2 for the oxidation of redox targets. Altogether, these results suggest that H2O2 signaling is mediated by localized redox relays whereby peroxiredoxins are oxidized to sulfenate and disulfide forms at H2O2 supply sites and these forms in turn oxidize the redox targets near these sites.

Mass isotopomer analysis of nucleosides isolated from RNA and DNA using GC/MS.

Nucleosides are biosynthesized from metabolites that are at key nodes of intermediary metabolism. Therefore, (13)C labeling patterns in nucleosides from ribonucleic acid (RNA) and deoxyribonucleic acid (DNA) in suitably designed isotopic tracer studies provide information on metabolic flux distributions of proliferating cells. Here, we present a gas chromatography (GC)-mass spectrometry (MS)-based approach that permits one to exploit that potential. In order to elucidate positional isotopomers of nucleosides from RNA and DNA, we screened the fragmentation spectra of their trimethylsilyl derivatives. We identified the molecular ion moieties retained in the respective fragment ions, focusing particularly on the carbon backbone. Nucleosides fragmented at the N-glycosidic bond provide nucleobase and/or ribose or 2'-deoxyribose fragment ions and fragments thereof. Nucleoside fragments composed of the nucleobase plus some carbons of the ribose ring were also observed. In total, we unequivocally assigned 31 fragments. The mass-isotopic distribution of the assigned fragments provides valuable information for later (13)C metabolic flux analysis as indicated by a labeling experiment applying [1-(13)C]glucose in a yeast culture.

Hydrogen peroxide metabolism and sensing in human erythrocytes: a validated kinetic model and reappraisal of the role of peroxiredoxin II.

Hydrogen peroxide (H2O2) metabolism in human erythrocytes has been thoroughly investigated, but unclear points persist. By integrating the available data into a mathematical model that accurately represents the current understanding and comparing computational predictions to observations we sought to (a) identify inconsistencies in present knowledge, (b) propose resolutions, and (c) examine their functional implications. The systematic confrontation of computational predictions with experimental observations of the responses of intact erythrocytes highlighted the following important discrepancy. The high rate constant (10(7)-10(8) M(-1) s(-1)) for H2O2 reduction determined for purified peroxiredoxin II (Prx2) and the high abundance of this protein indicate that under physiological conditions it consumes practically all the H2O2. However, this is inconsistent with extensive evidence that Prx2's contribution to H2O2 elimination is comparable to that of catalase. Models modified such that Prx2's effective peroxidase activity is just 10(5) M(-1) s(-1) agree near quantitatively with extensive experimental observations. This low effective activity is probably due to a strong but readily reversible inhibition of Prx2's peroxidatic activity in intact cells, implying that the main role of Prx2 in human erythrocytes is not to eliminate peroxide substrates. Simulations of the responses to physiological H2O2 stimuli highlight that a design combining abundant Prx2 with a low effective peroxidase activity spares NADPH while improving potential signaling properties of the Prx2/thioredoxin/thioredoxin reductase system.

Homeostasis of free cholesterol in the blood: a preliminary evaluation and modeling of its passive transport.

The rate of noncatalyzed transfer of cholesterol (Chol) among lipoproteins and cells in the blood is of fundamental importance as a baseline to assess the role of active transport mechanisms, but remains unknown. Here we address this gap by characterizing the associa-tion of the Chol analog, ergosta-5,7,9(11),22-tetraen-3ß-ol (DHE), with the lipoproteins VLDL, LDL, HDL2, and HDL3 Combining these results with data for the association of DHE with liposomes, we elaborated a kinetic model for the noncatalyzed exchange of free Chol among blood compartments. The computational results are in good agreement with experimental values. The small deviations are explained by the nonequilibrium distribution of unesterified Chol in vivo, due to esterification and entry of new unesterified Chol, and eventual effects introduced by incubations at low temperatures. The kinetic profile of the homeostasis of unesterified Chol in the blood predicted by the model developed in this work is in good agreement with the observations in vivo, highlighting the importance of passive processes.

Is Peroxiredoxin II's peroxidase activity strongly inhibited in human erythrocytes?

H2O2 elimination in human erythrocytes is mainly carried out by catalase (Cat), glutathione peroxidase (GPx1) and the more recently discovered peroxiredoxin 2 (Prx2). However, the contribution of Prx2 to H2O2 consumption is still unclear. Prx2's high reactivity with H2O2 (kPrx2=10×10(7) M(-1)s(-1), kCat =7×10(7) M(-1)s(-1), kGPx1 =4×10(7) M(-1)s(-1)) and high abundance ([Prx2]= 570µM, [Cat]= 32µM, [GPx1]= 1µM) suggest that under low H2O2 supply rates it should consume >99% of the H2O2. However, extensive evidence indicates that in intact erythrocytes Prx2 contributes no more than Cat to H2O2 consumption. In order for this to be attained, Prx2's effective rate constant with H2O2would have to be just ~10(5) M(-1)s(-1), much lower than that determined in multiple experiments with the purified proteins. Nevertheless, nearly all Prx2 is oxidized within 1min of exposing erythrocytes to a H2O2 bolus, which is inconsistent with an irreversible inhibition. A mathematical model of the H2O2 metabolism in human erythrocytes [Benfeitas et al. (2014) Free Radic. Biol. Med.] where Prx2 either has a low kPrx2 or is subject to a strong (>99%) but readily reversible inhibition achieves quantitative agreement with detailed experimental observations of the responses of the redox status of Prx2 in human erythrocytes and suggests functional advantages of this design (see companion abstract). By contrast, a variant where Prx2 is fully active with kPrx2=10(8) M(-1)s(-1) shows important qualitative discrepancies. Altogether, these results suggest that Prx2's peroxidase activity is strongly inhibited in human erythrocytes. We acknowledge fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) co-financed by FEDER through the COMPETE program and by FCT.

Is the Peroxiredoxin 2/Thioredoxin/Thioredoxin Reductase system in human erythrocytes designed for redox signaling?

In human erythrocytes H2O2 is mainly consumed by glutathione peroxidase, catalase and peroxiredoxin 2 (Prx2). Our previous analyses indicate that Prx2's peroxidase activity is subjected to a strong but quickly reversible inhibition (see companion abstract). If this activity is inhibited then the main role of Prx2 cannot be to eliminate H2O2. What functional advantages could then such an inhibition confer?We set up and validated a kinetic model of H2O2 metabolism human erythrocytes that shows quantitative agreement with extensive experimental observations. We then applied it to analyze the behavior of Prx2 and Trx under the H2O2 exposure dynamics that erythrocytes face in circulation. The significance of Prx2 inhibition was assessed by comparing the behavior of this model with that of an otherwise identical model lacking inhibition.Our analysis shows that Prx2 inhibition leads to 25-40% lower NADPH consumption under low to moderately high H2O2 supply (<0.8µM H2O2/s). Further, the inhibition extends the range where the concentrations of potential redox signaling readouts - H2O2, Prx2 sulfenic acid, Prx2 disulfide and Trx disulfide- show a proportional response to changes in H2O2 supply, covering practically the whole physiological range of the latter. This is desirable for analogic signal transduction and allows the Prx2/Trx/TrxR system to reliably transduce changes in H2O2 supply as changes in thiol oxidation. Finally, the inhibition allows other less abundant peroxiredoxins in the erythrocyte to be oxidized by H2O2 at physiological H2O2 supplies.Altogether, these results suggest that the postulated reversible inhibition of Prx2's peroxidase facilitates signal transduction by the peroxiredoxins and spares NADPH.We acknowledge: fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 co-financed by FEDER through the COMPETE program and by FCT (project PTDC/QUI-BIQ/119657/2010).

Integration between anticipatory blocking and redox signaling by the peroxiredoxin/thioredoxin/thioredoxin-reductase system.

Cells are occasionally exposed to high H2O2 concentrations, often preceding exposure to other electrophylic compounds. Both H2O2 and these compounds can irreversibly modify protein thiols, with deleterious consequences. Induction of enzymatic defenses against those agents is too slow to avoid significant damage. Cells may solve this conundrum by reversibly "blocking" the thiols once H2O2 concentrations begin to increase. We term this mechanism "anticipatory blocking" because it acts in anticipation of irreversible damage upon detection of early signs of stress. Here we examine the design requirements for the Peroxiredoxin/Thioredoxin/Thioredoxin-Reductase/Protein-Dithiol System (PTTRDS) to effectively integrate H2O2 signaling and anticipatory blocking of protein dithiols as disulfides, and we compared them to the designs found in cells. To that effect, we developed a minimal model of the PTTRDS, and we defined a set of quantitative performance criteria that embody the requirements for (a) efficient scavenging capacity, (b) low NADPH consumption, (c) effective signal propagation, and (d) effective anticipatory blocking. We then sought the design principles (relationships among rate constants and species concentrations) that warrant fulfillment of all these criteria. Experimental data indicates that the design of the PTTRDS in human erythrocytes fulfills these principles and thus accomplishes effective integration between anticipatory blocking, antioxidant protection and redox signaling. A more general analysis suggests that the same principles hold in a wide variety of cell types and organisms. We acknowledge grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) financed by FEDER through the "Programa Operacional Factores de Competitividade, COMPETE" and by national funds through "FCT, Fundação para a Ciência e a Tecnologia".

Organizing principles underlying microorganism's growth-robustness trade-off.

Growth Robustness Reciprocity (GRR) is an intriguing microbial manifestation: the impairment of microorganism's growth enhances their ability to resist acute stresses, and vice-versa. This is caused by regulatory interactions that determine higher expression of protection mechanisms in response to low growth rates. But because such regulatory mechanisms are species-specific, GRR must result from convergent evolution. Why does natural selection favor such an outcome? We used mathematical models of optimal cellular resource allocation to identify the general principles underlying GRR. Non-linear optimization allowed to predict allocation patterns of biosynthetic resources (ribosomes devoted to the synthesis of each cell component) that maximize growth. These models predict the down-regulation of stress defenses under high substrate availabilities and low stress levels. Under these conditions, stress tolerance ensues from growth-related damage dilution: the higher the substrate availability, the fastest the dilution of damaged proteins by newly synthesized proteins, the lower the accumulation of damaged components into the cell. In turn, under low substrate availability growth is too slow for effective damage dilution, and the expression of the defenses up to some optimal level then increases growth. As a consequence, slow-growing cells are pre-adapted to withstand acute stresses. Therefore, the observed negative correlation between growth and stress tolerance can be explained as a consequence of optimal resource allocation for maximal growth. We acknowledge fellowship SFRH/BPD/90065/2012 and grants PEst-C/SAU/LA0001/2013-2014 and FCOMP-01-0124-FEDER-020978 financed by FEDER through the "Programa Operacional Factores de Competitividade, COMPETE" and by national funds through "FCT, Fundação para a Ciência e a Tecnologia" (project PTDC/QUI-BIQ/119657/2010).